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Abstract 



We outline the most recent theory for the computation of the exponential 
growth rate of the number of configurations on a multi-dimensional grid. As 
an application we compute the monomer-dimer constant for the 2-dimensional 
grid to 8 decimal digits, agreeing with the heuristic computations of Baxter, 
and for the 3-dimensional grid with an error smaller than 1.35%. 
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1 Introduction 

The exponential growth rate h (with respect to the natural logarithm) of the number 
of configurations on a multi-dimensional grid arises in the theory of various phenom- 
ena [29, 10]. In physics is viewed as the entropy (per atom) of the corresponding 
"hard model"; in mathematics h is called the topological entropy [12]; and in infor- 
mation theory h (with respect to log2) is called the -multi- dimensional capacity [33]. 
The 1-dimensional case is easy, namely is equal to the spectral radius p{A) of 
a certain matrix A called the "transfer matrix". There are very few 2-dimensional 
models where the value of h is known in closed form [8, 21, 24, 25, 2]. In all other 
cases there are estimates based on: (a) asymptotic expansions, e.g., [27, 1, 15]; (b) 
Monte-Carlo methods, e.g., [18, 3]; (c) bounds, e.g., [17, 7, 26, 6, 9, 28]. In what 
follows we give a complete up-to-date theory of the computation of h by using lower 
and upper bounds. It refines the techniques described in [13] by using an automor- 
phism subgroup of a given graph. A fundamental problem in lattice statistics is the 
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monomer-dimer problem (see [22]). As a demonstration of our techniques, we com- 
pute the topological entropy of the monomer-dimer covers of the 2-dimensional grid 
/i2 = .66279897 (which agrees with the heuristic estimation e^^ = 1.940215351 due 
to Baxter [1]) and of the 3-dimensional grid .7653 < /is < .7862. These numerical 
results are much better than previously known ones. 

Consider the grid Z"^ in d-dimensional space M"^. At each point of the grid we 
place an element of a set of n kinds of colors (atoms) denoted by (n) := {1, . . . , n}. 
Certain restrictions may be imposed on the colorings. For example, the restrictions 
of the hard model are specified by a directed d-graph F := (Fi,...,F(^) called a 
nearest neighbor digraph, with F^ C (n) x (n), in the sense that two atoms of kinds 
p and q are allowed to occupy respectively the adjacent grid points i = (?i, . . . ,1^) 
and i -|- e^. (where := {Sik, ■ ■ ■ , ^dk)) only if (p, q) G F/;. We call such a placement 
a T -configuration or T-cover. This general model is anisotropic, since the F^ can 
be distinct. A digraph F^ is called symmetric when {p.,q) G F^ 44> (QiP) G F^. We 
call F a symmetric isotropic nearest neighbor digraph when Fi = • • • = F(i = A, and 
A is symmetric. Let m = (mi, . . . , m^) £ N*^, where N := {1, 2, . . .}, and consider 
the box (m) := (mi) x ••• x {md) of dimensions mi, . . . ,md. Let Ty(m) be the 
set of all F-configurations of |m|pi- := mi • • • m^ atoms in the box (m). It is easy 
to show that the sequence log #VF(m)^gj^d is subadditive in each coordinate, i.e., 
log #VF(m-|-pefc) < log#VK(m)-|-log #W{m+{p—mk)ek) for all m G N"', p €N and 
k £ (d). From this it can be shown that the following limit exists and is non-negative 
or equal to — oo (we use m ^ oo as an abbreviation of mi, ... , ma — > oo): 

, = ,(r):= l.m !55#2>^^), (1.1) 

m^oo |m|pr 

and each m G N*^ satisfies 

,.<'°g*f'-). (1.2) 

|m|pi. 

The limit h{T) is the exponential growth rate of i^W{m) per atom, also called 
entropy or Shannon capacity. It follows from Konig's Infinity Lemma that h = 
logO = — oo if and only if there are no F-covers of Z'^. The case d = 1 is well 
understood: h = log p{A), where A is the incidence matrix for the digraph Fi; 
there exist F-covers if and only if Fi has a directed cycle, and in that case h is 
also the exponential growth rate per atom of the number of periodic F-covers of Z 
[12]. A periodic F-cover of Z'^ with period m (i.e., a F-cover cf) = {(t)i)i^z'i satisfying 
<^i+mfeefc = 0i all i G Z"^ and k G (d)) is equivalent to a F-cover of the torus 
T(m) := (Z/miZ) x ••• x (Z/m^Z). For d > 2, the question whether there exist 
F-covers is undecidable and h is not computable in general [4, 20] (we say that a 
quantity Q is computable when given e > 0, we can find in a finite number of steps, 
depending on e, a rational number r satisfying |Q — r| < e). Equivalently, there is a 
d-digraph F for which there are F-covers of Z"^ but none is periodic. Hence there are 
no nontrivial lower bounds for h in this case. A fundamental result in [12] asserts 
that if at least d — 1 digraphs out of Fi, . . . , F^ are symmetric, then the exponential 
growth rate per atom of the number of periodic configurations is equal to h and h 
is computable, i.e., we have lower bounds on h that converge to h. For d = 2,3 this 
will also follow from our results in Section 3. In particular these results hold for a 
symmetric isotropic nearest neighbor digraph. 
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We mention briefly the topological entropy. Let Wtop(ni) be the set of all distinct 
restrictions of F-covers of Z*^ to the box (m). log 7^Wtop(in) is also subadditive, and 
the topological entropy of T is defined by 

log#Tytop(m) 
'itop(r) := lim ■ — ■ . 

Since Wtop(m) C W{m), we have /itop(r) < ^(r); a result in [12] asserts that 
equality holds. 

We now elaborate our results. Fix m' := (mi, . . . ,mrf_i) G N'^""'^ and let V := 
(Fi, . . . ,Frf_i). Let Orf(m') be the transfer digraph between F'-covers of (m') with 
respect to F^. That is, the vertex set of Orf(m') is the set of F'-covers of (m'), and 
vertices u,v satisfy {u,v) G Qrf(m') if and only if [u,v] £ W{m',2), where [u,v] 
is the configuration consisting of u,v occupying the levels = 1,2 of ((m',2)), 
respectively. We show that h < ^°^^^,|^'" (by definition, the spectral radius 
of a digraph is the spectral radius of its incidence matrix). When Fi, . . . ,Frf_i are 
symmetric, this upper bound can be improved as follows. Let Gd(m') be the induced 
subdigraph of J7(;(m') whose vertices are the periodic F'-covers of (m') with period 
m'. Then we show [13] 

hfr\ ^ logp(ed(m')) r r + • n Q^ 

n(i ) < j — , mi, mrf_i even, 1 1, 1 symmetric. (1-3) 

|m Ipr 

Furthermore, for Fi, . . . , T^-i symmetric, we give various lower bounds on h in terms 
of log p(G(i(m')) for various values of m'. For example, for d = 2 we show that when 
Fi is symmetric, h > l°gp(e2(p+2g»-logp(e2(2g)) ^11 p G N and G Z+ := N U {0}. 
See [13] for slightly different lower bounds on h, which do not use periodicity. 

All of these upper and lower bound converge to the true entropy when m' — > oo. 

We can enhance the efficiency of computing the spectral radius /o(A) of a digraph 
A C X as follows. To compute p(A) one needs to compute the spectral radius 
of its 0-1 N X N incidence matrix A. Suppose that Q C Sn is an automorphism 
subgroup of A. Let O = {N)/Q be the orbit space under the action of Q and set 
M = #0. Let A' C O X O be the multidigraph induced by A and Q. That is, for 
a, f3 £ O, the multiplicity of the edge (q, (3) of A' is aa,/3 = J2je(3 i £ a. 

We show that p(A) is also the spectral radius of the M x M nonnegative integer 
matrix A. If M <C A^, then the computation of p{A) may be feasible on a desktop 
computer whereas the computation of p{A) may be infeasible on a supercomputer. 

We show that that the automorphism group of Grf(m') contains a subgroup 
isomorphic to the group of translations of T(m'). If Fi = • • • = F^^i = A and A is 
symmetric, then the automorphism group of ©^(m') contains a subgroup isomorphic 
to the group of rigid motions of T(m') (motions preserving the distance on T(m'), 
i.e., translations, reflections and coordinate transpositions for equal dimensions). 
For example, T{m) has m translations and 2m rigid motions if m > 2. 

We now discuss the monomer-dimer covers of Z'^, see [10]. A diraer is a domino 
consisting of two neighboring atoms occupying the places i, i + e/j G Z'^. A monomer 
is a single atom occupying the place i G Z'^. A monomer-dimer cover, respectively 
dimer cover, of Z"^ is a partition of Z*^ into monomers and dimers, respectively 
dimers. We denote by hd and the entropies of the monomer-dimer and dimer 
covers, respectively. It is fairly easy to compute the values hi = log and hi = 0. 
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The big breakthrough in the sixties was a close formula for /i2 in [8, 21]. The exact 
values of for d>2 and /i^ for d> 3 are unknown. 

A seminal contribution to the study of upper and lower bounds and estimates 
for hd and hd was given in [16, 17, 18, 19]. In particular, it was shown in [16] that 
for p G [0,1], there exists the entropy Xd{p) of the monomer-dimer covers of Z'^, 
where p is the "density" of dimers, i.e., the number of dimers in the cover divided 
by one half of the volume. The entropy Xd{p) is a continuous concave function of p 
and Ad(l) = hd- It is shown here that hd = maXpg[o^i] Xd{p)- It was pointed out by 
Kingman, see [17], that the van der Waerden conjecture for permanents of doubly- 
stochastic matrices gives a lower bound on hd- The improved lower bound for the 
permanents of 0-1 matrices [31] gives the currently best lower bound /13 > 0.440075. 
A recent breakthrough [7] gives the upper bound /13 < 0.463107, improved in [26] 
to hs < 0.457547. 

It is shown in [13] that the dimer covers can be encoded as A-covers for an 
appropriate d-digraph A = (Ai, . . . , Ad), where all digraphs are on the set of vertices 
(2d). We show that the monomer-dimer covers can be similarly encoded as A- 
configurations for an appropriate d-digraph A = (Ai,...,Arf), where all digraphs 
are on the set of vertices (2d -|- 1) . Unfortunately, in these encodings the digraphs 
Ffc, Ffc are not symmetric, so (1.3) and the lower bounds do not apply directly. One of 
the purposes of this paper is to show that the entropies hd and hd nevertheless obey 
upper and lower bounds converging to the true entropies, similar to (1.3) and the 
lower bounds discussed above for the symmetric isotropic nearest neighbor digraph. 
The bounds for hd are stated in terms of the spectral radii of certain multidigraphs 
0d(m') whose automorphism group has a subgroup isomorphic to the the group of 
rigid motions of T(m'). This fact enables us to compute the values of /i2 and /13 
with good precision. We also show that Xd{p) can be bounded below by using the 
generalized van der Waerden conjecture (Tverberg's conjecture), proved by the first 
author in [11]. For d = 2, 3, this lower bound is better than those of [5] and [19] 
except for very high p. Our lower bound for Xd{p) yields in particular a lower bound 
for hd- For d = 2 this lower bound is somewhat weaker than the one obtained from 
the numerical computations of p(Gc;(m')), but for d = 3 the situation is reversed. 

See [14] for a general theory of monomer-dimer covers of an arbitrary graph. 
Finally it is worth mentioning the theoretical work [23] , which shows that the general 
monomer-dimer problem in arbitrary planar graphs is computationally intractable. 

The contents of the paper is as follows. In Section 2 we discuss the the general 
theory of Z"^ subshifts of finite type (SOFT). In Section 3 we prove the main in- 
equalities of the entropy of Z'^-SOFT with d — 1 symmetric digraphs Fi, . . . ,Td-i- 
In Section 4 we recall the main features of the entropy of the monomer-dimer and 
dimer covers. In Section 5 we give lower bounds for the entropy of the monomer- 
dimer covers with a fixed dimer density using the lower bounds on permanents. In 
Section 6 we show that there exist analogs of the upper and lower bounds discussed 
in Section 3 that apply to the monomer-dimer and dimer entropy. In Section 7 we 
discuss using automorphism subgroups to reduce the computations. In Section 8 we 
give numerical upper and lower bounds for /i2, /i2i h^, /13, and compare graphically 
our lower bounds for A2(p) and X-^lp) with the known lower bounds and estimates. 
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2 SOFT and NNSOFT 



Let (n)^"^ be the set of ah colorings (j) : 'L'^ ^ (n) of with colors from {n) = 
{l,...,n}. Given a d-digraph T = {Ti,...,Td) on (n) x (n), let T^'* C (n)^'* be 
the set of all F-covers, namely colorings (j) = (0m)meZ'* ("')^ such that for each 
1^1/ and k G (d) , the restriction of (f) to the line through i in the direction of , 
i.e., (</>i+jefe)jeZ5 is a bi-infinite walk on F^. In ergodic theory, F^'* is called a nearest 
neighbor subshift of finite type (NNSOFT). Note that for an NNSOFT T^'' and for 
m G N'^, VF(m) is the set of all configurations G (n)^""^ such that i, i + G (m) 
imply (V'i,V'i+eJ G ^k- 

A general SOFT can be described as follows. Let M G N'^ and a nonempty 
subset V C (n)^'^^ be given. Every element a G "P is viewed as an allowed coloring 
(configuration) of the box (M) in n colors. For i G Z'^, we define the shifted coloring 
Ti{a) of a G "P as the coloring of the shifted box (M) + i that gives to the point x + i 
the same color that a gives to x G (M). We denote by Ti('P) the set {ri(a) : a G V}, 
and regard it as the set of allowed colorings of (M) + i. A coloring (f> G (n)^ is 
called a V -state if for each i G Z*^ the restriction of to (M) + i is in T{(V). We 
denote by (n)^ (V) the set of all P-states. In ergodic theory the set (n)^ (V) is 
called a subshift of finite type (SOFT) [30]. 

Each NNSOFT V'^" is a special kind of SOFT obtained by letting M = (2, . . . , 2) 
and V the set of all colorings ijj G (n)^'^^ such that i, i+e^ G (M) imply (-(/'i , -i/'i+efe ) £ 
Ffc. Conversely [12], each SOFT {n)'^\v) can be encoded as an NNSOFT F^^ 
where F = (Fi, . . . ,Frf) are defined as follows. Take = #P and use a bijection 
between V and (A^). The digraph F^ C (A^) x (A^) is defined so that for a,b € V we 
have (a, 6) G F^ if and only if there is a configuration (p G (n)^^^'^'''^ such that the 
restriction of cp to (M) is a and the restriction of cp to (M) + is Te^{b). Because 
of this equivalence, we will be dealing here with NNSOFT only. 

In the sequel we will be taking lim sup and lim inf of real multisequences (am)mGN'* 
as m ^ oo. In order to be clear, we define these here and prove that they are limits 
of subsequences. We also define the limit of real multisequence in terms of lim sup 
and lim inf, which is equivalent to other definitions in the literature. 

Definition 2.1 Let (am)mGN<' ^ multisequence of real numbers. Then 

(a) limsupm^j^ Om is defined as the supremum (possibly ±ooJ of all numbers 
of the form limsupg^g^ am,, where (mg)ggM is a sequence in N'^ satisfying 
limq^oo = OO, i.e., limq^oo(mg)j = oo for each i G (d). We define 
liminfin-»oo fflm similarly. 

(b) limm_^oo flm = OL means limsupj^^oo = liminfjn-»oo Om = «• 

Proposition 2.2 //limsupj^_^oQ Om = ot, then there exists a sequence (nq)ggN C 
N*^ satisfying linig^oo = oo such that the sequence (ang)qeN has a limit and 
liniq^oo = Oi. Similarly for lim inf. 

Proof. Since the lim sup of each real sequence is the limit of a subsequence, 
we may assume that we have a sequence of convergent subsequences {om* } satisfying 
limq__»oo Om* = cti ^i^d limq_»oo = oo for each i G N, and limj^oo = Note 
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that ai < a for all i by definition of the supremum. If = a for some i, we are 
done. In particular, if a = —oo, then a, = — oo = a for each i and we are done. 
Therefore we may assume that a G M U {00} and that ai < ctj+i for all i. 

Assume first that a S M. Then for each i E N there exists a q{i) £ N such 
that m^~l^,-,-, > 2m%.s and la^i — Ojl < ^- Then we can take rij = m*,.N and 

the result follows. Similarly, if a = 00, then for each i E N there exists a q{i) E N 
such that m^"^.^j^^ > 2m* and a^i ^ ^ > ai — 1, and again we can take rij = ^ 

Let Typer(ni) C F^'' be the set of periodic F-covers with period m. Then 

nper(r) := limsup j — p (2.1) 

m^oo |m|pr 

is called the periodic entropy of F^'' . Clearly /iper(r) < h{T). 



3 Main Inequalities for Symmetric NNSOFT 

For d > 2, consider m = (mi, . . . , m^) E N'^ and m~ := {m2,---,md). Let 
Wpgj, |ij(m) be the set of F-configurations in the box (m) that correspond to F- 
covers of T{mi) x (m~), i.e., that can be extended periodically in the direction of 
ei with period mi into F-covers of Z x (m~). We can view these configurations 
as F-configurations in the box (m~), where F = (F2, . . . ,Frf), for each j the vertex 
set of Tj is the set F^p^^^, of closed walks a = (ai, . . . , a^i , oi) of length nii on Fi, 

and where (a, b) E Fj if and only if (oj, 6j) E Tj for i = 1, . . . , mi. For this reason, 
the following limit exists and is equal to the entropy /i(F2, . . . , F^) of the NNSOFT 

n(mi,F) := lim ■ — _ , mi E N (3.1) 

m~^oo |in Ipi- 

We define W~{m~) as the set of (F2, . . . ,Frf)-covers of the box (m^). In the degen- 
erate case mi = 0, we define Wper^{i}.(0, m~) to be simply W~{m~) and (F2, . . . ,Frf) 
to be simply (F2, . . . , F^). Then (3.1) is also valid for mi = 0, where we understand 
MO,F) to be /i(F2,...,Fd). 

Theorem 3.1 Consider the NNSOFT T^" for d > 2. Let h{T) and h{r,T) be 
defined by (1-1) and (3.1), respectively. Assume that Fi is symmetric. Then for all 
p, r E N and q E Z+, 

M2r£) ^ ^^^^ ^ % + 2,,F)-M2g,F)^ ^3^^^ 
2r p 

Proof. Fix m^ = (m2,...,mrf) E N'^^^ and let Oi(m~) be the following 
transfer digraph on the vertex set W~ {in~), analogous to the transfer digraph 
rirf(m') described in Section 1. Vertices u,v satisfy {u,v) E rii(m~) if and only 
if [u,v] E W{2,m~), where [u,v] is the configuration consisting of u,v occupying 
the levels xi = 1,2 of ((2,m^)), respectively. Let N = (m^) and let C(m~) 

be the x 0-1 incidence matrix of r2i(m~), with spectral radius p{C{in^). As 
a nonnegative matrix, C(m~) satisfies (see e.g., [13]) 

logp(G(m )) = Inn ; , 

fc— >oo k 
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where 1 = (1, . . . , l)"*". Since l~^C{m is the number of walks of length k on 
r2i(m~), which correspond to F-covers of ((/c,m^)), we obtain 

logp(C(m')) ^ log#l^(fc,m-) 
|m~|pi. fc^oo fc|m~|pi. 

Now send m2, . . . , ma to oo, and observe that by (1.1) and (1.2), the right-hand side 
of (3.3) converges to h{T) and is an upper bound on it for each m~. Thus we obtain 
[12] 

^°gf(^^,"^"^) > Mr), m- G (3.4) 

|m Ipr 

hm i^^^M^ = Mr). (3.5) 



m -^oo in 



Ipr 



Next, we observe that 

trC(m")'? = #T^per,{i}(<7,m"), g G Z+, (3.6) 

where C{m~)^ is the N x N identity matrix. Recall that the trace of C{m~)'^ is 
given by 

N 

trC(m-)^ = ^Af, 
1=1 

where Ai,...,A7v be the eigenvalues of C(m~). Since C(m~) is a nonnegative 
matrix, the Perron- Frobenius theorem yields that its spectral radius p{C{m~)) := 
maxjg^^^ |Aj| is one of the Aj. Since by assumption Fi is symmetric, r2i(m~) and 
hence C(m~) are symmetric. Therefore Ai, . . . , Atv are real, and hence tr C(m~)^'" > 
p(C(m~))^'' for each r G N. Taking logarithms and using (3.6), we obtain 

log#^pcr.{i}(2r,m') ^ logp(C(m-)) ^ ^ ^ 
2r|m-|pi. ~ |m-|pr 

Sending m2, . . . ,mrf to oo in (3.7) and using (3.1) and (3.5), we deduce the upper 
bound on /i(F) in (3.2). 

To prove the lower bound in (3.2), we note that 

trC(m-)P+2'? = J] Af+29 < \^^\'"^^' = Yl l^^l"^*'' 

i i i 

< J^p(C(m"))PA,2'' = p(C(m^))PtrC7(m")29 

i 

and thus by (3.6) 

logp(C(m-)) ^ log # Wpcr.ji} (p + 2g, m" ) - log # Wper,{i} (2g, ) 
|m-|pr ~ p\m~\pr 
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Sending m to oo and using (3.5) and (3.1) (recall that the latter holds for mi G Z+), 
we deduce the lower bound in (3.2). □ 

When d = 2, h{mi,T) is the entropy of the NNSOFT ff (recall that is 
simply r2 when mi = 0). Since this is a 1-dimensional NNSOFT, that entropy is 
equal to logp{T2)- We denote p{T2) by 62{rni), and obtain the following corollary 
to Theorem 3.1. 

Corollary 3.2 Let d = 2 and assume that Ti is symmetric. Then for allp, r G N 
and q G Z_|_, 

log^2(2r) ^ ^^^^ ^ logg2(p + 2(z)-logg2(2g) ^ ^3 
2r ~ ~ p ' 

where 02 is defined above. 

In (3.9) take q = and p = 2r, and send r to oo. Clearly the upper and lower 
bounds then converge to h{T). Hence h{T) is computable [12]. For completeness of 
the exposition we reproduce a short proof of (1.3) for any d > 2 given in [13]. We 
use the following straightforward lemma. 

Lemma 3.3 Let T = (Fi, . . . .F^) and m £ N"', put F' = (Fi, . . . ,Td-i) and 
m' = {mi, . . . ,md-i), and let 0d(m'J be the transfer digraph between T' -covers of 
T(m') with respect to F^. Let F2, . . . , F^ 6e defined as in the beginning of this section, 
put T' = (F2, . . . ,Td~i) and m = (m2, . . . ,r?irf_i), and let ©^(m) be the transfer 
digraph between T' -covers ofT{m) with respect to F^. Then 0^(m') and ©^(m) are 
isomorphic, and in particular /9(0^(m')) = /9(0rf(m)). 

Proof. We use the following bijection between the vertices u of 0(i(m') and the 
vertices u of 0(i(m). Given u = (0i)i(=(m')5 '^^ have 

(0i,(/)i+ej G Ffc, k=l,...,d-l, (3.10) 

where the addition i + is understood modulo m^, i.e., + 1 is 1. Then the 
corresponding u is defined to be n = ('/'j)jG(m)5 where (pj = {4>{q,}))q=i- We note that 
(pj is indeed a Fi-cover of T{mi) and thus a vertex of F2, . . . ,T(i-i by (3.10) with 
i = (5, j) and k = 1. In order to show that n is a F '-cover of T(m) and thus a vertex 
of 0d(m), we need to show that ((/>j, '/'j+e'^.) ^ for k = 2, . . . ,d — 1. This means 
showing that (</'(gj), (/'(gj+e'^.)) G F^. for /c = 2, . . . , d — 1 and q = 1, . . . , mi, which 
follows in turn from (3.10) with i = ((/,j). It is easy to see that the correspondence 
u^u can be inverted. It remains to show that {u,v) G 0(i(m') 44> {u,v) G 0rf(m). 
We prove only the =^ part. Let u = {4>i)\^{m') ai^d v = (^i)ig(m') be F'-covers of 
T(m'). The assumption {u,v) G 0rf(m') means that G F^ for all i G (m'). 

Applying this with i = (g, j), q = 1, . . . , mi and j G (m) shows that ((/)j, ^j) G F^ for 
all j G (m), which means in turn that (u,v) G 0(i(m). □ 



Theorem 3.4 Letd>2 and consider the NNSOFT , where F = (Fi, . . . , F^). 
For ui' = {mi, . . . ,mrf_i) G N'^""'^ and T' = (Fi, . . . ,Td-i), let 0rf(m') be the trans- 
fer digraph between V -covers o/T(m') with respect to F^. Assume that Fi, . . . , F^^i 
are symmetric and mi, . . . , m^^-i are even. Then 

< iogp(e^. 

jm'lpr 
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Proof. The proof is by induction on d. For d = 2 the result is equivalent to the 
upper bound in (3.9). For the induction step, observe that the upper bound of (3.2) 
with r = mi/2 yields h{T) < h{mi,T)/mi. Recall that h{mi,T) is the entropy of 
the NNSOFT " , where F = (Fg, . . . , Fd) is as in Lemma 3.3. Since F2, . . . , F^.i 
are symmetric, so are F2, . . . ,Fd_i, and therefore the induction hypothesis applied 
to F^ gives /i(mi, F) < log p(0rf(m)) /|m|pr, where m and are as in the lemma. 
Finally, an application of the lemma completes the proof. □ 



Corollary 3.5 Let F = (Fi,F2,F3) and assume that Fi and F2 are symmetric. 
For (mi, 7712) G N^, let ©3(^1,7712) be the transfer digraph between (Ti,T 2) -covers of 
r(mi, m2) with respect to F3, and let 03(mi, 7712) be its spectral radius. Let ^3(0, 7712) 
be the spectral radius of the transfer digraph between T2-covers ofT{m2) with respect 
to F3. Let ^3(7711,0) be the spectral radius of the transfer digraph between Ti-covers 
ofT{mi) with respect to T^. Then for all r,t,p,u,v £ N andq,s£Z^ we have 

^oM^ > ftiD 

Art - ^ I 

> log 9^{p + 2q,u + 2s) - log 03 (P + 2g, 2s) _ logg3(2g, 27;) 
~ up 2vp 

Proof. The upper bound in (3.11) follows directly from Theorem 3.4 for d = 3. 
To show the lower bound we use the lower bound in (3.2), which is valid since Fi is 
symmetric, and gives 

/.(r..r..r3)>^"' + ^''-'^"^""^'>-^P''-'^"^""'». (3.12) 

p 

For each a G N we have h[a, (Fi,F2,F3)) = h{T2,Tz), where F2,F3 are digraphs on 
the vertex set Tl as in the beginning of this section. Since F2 is symmetric, so is 
F2, and so we can apply the lower bound of Corollary 3.2 to (F2,F3) to obtain 

h(r2, r3) > l"g^3(a,n + 2.)-logg3(a,2.)^ 

u 

where d3{a,b) = ^(©3(6)) = /9(©3(a, 6)) by Lemma 3.3. Inequality (3.13) is also 
valid for s = 0, since we defined ^3(0,0) to be the spectral radius of F3, exactly as 
in Corollary 3.2 for the degenerate case. Using (3.13) for a = p + 2q gives 

h( ^0 (r r ^^ ^ logO-jip + 2q,u + 2s) - logOajp + 2q,2.s) 

+ 2g, (Fi,F2,F3)) > . (3.14) 

u 

Apply the upper bound of Corollary 3.2 to (F2,F3) to obtain 

Mf2,f3)< ^"^'f'^^^ (3.15) 
27; 

Inequality (3.15) is also valid for a = 0, since in that case (F2,F3) = (F2,F3), by 
Theorem 3.4 applied to (F2,F3), and by the definition of ^3(0, 27;). Using (3.15) for 
a = 2q gives 

77^0 (r r ^^ ^ log 6*3 (2g, 27;) 

^(2g, (Fi,F2,F3)) < . (3.16) 

27; 

Finally, substitution of (3.14) and (3.16) in (3.12) yields the lower bound of (3.11). □ 
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(a) (b) (c) 

Figure 1: (a) T-coloring of (m) = ((3,3)); (b) Corresponding monomer-dimer cover 
of T; (c) Corresponding monomer-dimer cover of S 

4 Dimer and Monomer-Dimer Covers of 

As in [13], the set of monomer-dimer covers, respectively dimer covers, of Z'^ is an 
NNSOFT , respectively , where F and F are defined as follows. We encode a 
monomer-dimer cover of 'Z'^ as a coloring of 'Z'^ with the 2d+ 1 colors 1, . . . ,2d + 1: 
a dimer in the direction of e,fc occupying the adjacent points i, i + is encoded by 
the color /c at i and the color k + d at i + e^; a monomer at i is encoded by the color 
2d + 1 at i. This imposes restrictions on the coloring, which are expressed by the 
d-digraph F = (Fi, . . . , F^) on the set of vertices {2d + 1), where 

• {k, q) e Tk <^ q = k + d] 

• for j / k, (j, q) £Tk <^ q ^ k + d. 

It is easy to check that this gives a bijection between the monomer-dimer covers of 
Zi'^ and F^ . Similarly, if F = (Fi, . . . , F^^) is obtained from F by removing the vertex 
2d + 1, then there is a bijection between the dimer covers of Z*^ and F^''. 

The disadvantage of these encodings is that F^ and F^ are not symmetric, so 
we cannot apply the results of Section 3 directly. However, as pointed out in [7] for 
the dimer problem, there is a hidden symmetry, which enables us to obtain results 
analogous to those of Section 3. 

Recall that Vl^(m) denotes the set of F-colorings of (m) C N'^. Consider a F- 
coloring (j) £ W{m) with the F defined above. Certain points i on the boundary 
of (m) can receive colors indicating that i is one half of a dimer whose other half 
is outside (m). Therefore (p corresponds to a monomer-dimer cover of a "box with 
protrusions" T satisfying (m) C T C (m + 21) — 1, where 1 := (1,...,1) G N'^, such 
that each monomer in the cover is contained in (m) and each dimer in the cover has a 
nonempty intersection with (m). We translate T by 1 to move it into N*^, and thus (f) 
corresponds to a monomer-dimer cover of a set S satisfying (m) + 1 C 5" C (m + 21) 
such that each monomer in the cover is contained in (m) + 1 and each dimer in the 
cover has a nonempty intersection with (m) + 1. Conversely, each monomer-dimer 
cover of such a set S satisfying these conditions corresponds to a F-coloring of (m) . 
This is illustrated in Figure 1. 
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Similarly, W(m) denotes the set of F-colorings of (m), and there is a bijection 
between VF(m) and the set of dimer covers of a set S satisfying (m) + 1 C 5 C 
(m+21) such that each dimer in the cover has a nonempty intersection with (m) + l. 

Let Wper(m), respectively Wper(m), denote the set of F-colorings, respectively 
F-colorings, of (m) that can be extended periodically to F-colorings, respectively F- 
colorings, of Z'^ with period m. It corresponds to the set of monomer-dimer covers, 
respectively dimer covers, of T(m) and satisfies VFpcr(ni) C W{m), M^pcr(in) C 
Wim). _ 

Finally, let Wo(m)) respectively Wo{in), be the set of F-colorings of (m) for 
which S defined above is equal to (m) + 1, i.e., each dimer in the corresponding 
cover of S is contained in (m). To emphasize the fact that the dimers do not 
protrude out of (m), we refer to these covers as tilings. We have VFo(m) C Wpcr{i^), 
Wo{m) C 14/pcr(m). We can see that #VF(m) < i^Wo{m + 21), because we can 
extend the monomer-dimer cover of S into a member of Wo((m + 21)) by tiling 
(m + 21) \ S with monomers. 

From the discussion above we have 

#Wo(m) < #WVer(m) < #W{m) < #Wo{ui + 21) (4.1) 
#Wo{m) < #WpUm) < #W{in) (4.2) 
#Wo{m) < #Wo{ui) (4.3) 

#Wper{m) < #VFper(m) (4.4) 

#W{m) < #W{m). (4.5) 

Recall that the d-dimensional monomer-dimer entropy is defined by 

log#VF(m) 



hd := lim 



From (4.1) we obtain 



. ^ log #1^0 (m) . log#t^o(m + 21) 

lim mr ; — ; = lim mf ■ 



|m|pr m-+oo |m + 21|pr 



. ^ log#T^o(m + 21) ^ ^ log#VFo(m) 
hm mf j — j > h(i> hm sup j — j , 



This and one more application of (4.1) give 

\og#W{m) _ ^. log#l^per(m) 



hd '■= lim ; — ; = lim 



m-ioo Imlpr m^oo |in|pi. 



Inn '°'^f':°'-'. (4.6) 



Similarly, the d-dimensional dimer entropy is defined by 

log #W{m) 



hd := lim 

m— >oo 
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It is known to satisfy 

ha-= lim l°g#^("^) = lirn lQg#^Pcr(m) 

■ m^oo |m|pr „^oo,^eN l^lpr 



lim '"'^fr-'-l (4^7) 



The proof of (4.7) is more involved, and follows from the results proved in [16], as 



we show now. For m S N and s E 



0. 



|m|pr 



n Z, let Wo(m, s) be the subset of 



2 

VFo(m) consisting of the monomer-dimer tilings of (m) that have exactly s dimers. 
As pointed out in [16], WQ{m,s) 7^ by induction on d. It is shown in [16] that 
there exists a function Xd{-) '■ [0, 1] — > such that for all sequences (mq)qgN and 
{sg)g(zn satisfying 



Sg G 



0, '"'^'p'' 



2 

the following equality holds 



2s 

nZ, lim m„ = cx), lim- ^=pG[0, 1], (4.8) 

(?— >oo g— >oo jiriqlpr 



log#VFo(mg,Sg 



lim "\ = Up)- (4.9) 

q^oo Imqlpr 

Furthermore, the function Xd{p) is a continuous concave function of p on [0, 1]. We 
call Xd{p) the monomer-dimer entropy with dimer density p. 

Theorem 4.1 Let W{m), Wper(m), tVo(m) be defined as above. Then (J^.l) 
and the following equalities hold 

Xd{0) = (4.10) 
Xd{l) = hd (4.11) 
max Xdip) = hd. (4.12) 

PG[0,1] 

Proof. The proof of (4.10) is easy: pick any sequence niq satisfying linig^oo = 
00, and take Sg = for all q. Then conditions (4.8) hold for p = 0, and consequently 
(4.9) holds. But #lVo(mg,0) = 1, since there is only one way to cover a box with 
monomers, and (4.10) follows. 

We prove (4.7) and (4.11) together. Pick a sequence (mq)ggN ^1 N'^ such that 
the I 

Hiqlpr are even and limg_>oo 

nig = 00, and take Sg = Then conditions 

(4.8) hold for p = 1, and consequently (4.9) holds. But Wo{mq,Sq) = Wo{mq), and 
therefore 

lim »M = A.(1). (4.13) 

m^oo.l^eN |m|pr 

In view of (4.2) and (4.13) we obtain 

log#W>i) log#T^pcr(m) 
lim j — j > imisup j — f- 



> liminf l5i#i5WW> (4,14) 
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For m > (2, . . . , 2) G N*^ let a(m) := 2| 

'^Ipr Xyi=i rn~ ^"^^ surface area of (m), 



rrii 

and let 



w{m) := ^ #Wo(m,s) 



ix'(m) := max ^Wo{in,s) 



|m|pi— a(m) |m|prj 



s(m) := argmax #T¥o(m,s). 

.g[ l"'lp7'''"'\ ^]nz 

In words, i(;(m) is the sum of i^W(){m,s) where s ranges over those numbers of 
dimers that are sufficient to cover the interior of (m), i.e., the elements of (m) not 
on its boundary; s(m) is the largest summand in that sum; and ^(m) is a number 
of dimers achieving the maximum. 

Clearly ci;(m) < w{m) < °^"^^"'"^ ci;(m), and therefore 

\ogu{m) logu;(m) 
hmsup — . — . =limsup — . — . . (4.15) 

m— >oo l-'^lpr m— »oo l-'^lpr 

By Proposition 2.2 there exists a sequence (ng)ggN ^ satisfying 

r r log'^K) 1- log5(m) 
lim Uq = CO, hm — j — j — — = limsup — j — j . (4.16) 



q—*OD q^oo n, 



gipr 



Let tq := s{nq) for each g S N, and so ^Wo{iiq,tq) = uj{nq). Clearly limg^oo jirf~ ~ 
1, and so conditions (4.8) hold for nq,tq with p = 1, and consequently (4.9) holds 
for them. Hence by (4.16) 

limsupi^f^ = A,(l). (4.17) 

Next we assert that ^W{ni) < w{m + 21). Indeed, each cover in VF(m) can be 
shifted by 1 and extended by monomers to a tiling in IVo(iii + 21, s) for one of the 
s appearing in the sum w{m + 21). Therefore by (4.15) and (4.17) 

log#M^(m) log ii;(m + 21) 
lim j — j < lim sup j — j 



log,w(m) , logwfm) , , . 

= hmsup ^ ^ = limsup , , ^ = Ad(l). (4.18) 

m— >oo |m|pr m— too |m|pr 

Inequalities (4.14) and (4.18) combined, along with (4.2), complete the proof of (4.7) 
and (4.11). 

We now prove (4.12). As VFo(m, s) C Wo{ni), it follows that Xdip) < for 
all p S [0, 1]. To complete the proof, we exhibit a p* G [0, 1] satisfying the reverse 
inequality. For each m G N'^, let 

ti;(m) := max #Wo{m,s) 

se[o,^]nz 

s(m) := argmax ^Wo{TLn,s) 

se[o,l^]nz 

p{m) := e [0,1], 
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so that Lo{m) = ^Wo{m, s(m)). 

Observe t] 
fore, by (4.6), 



Observe that #Wo(m) = E^g^ jHprj^^ #VFo(m, s) < ■l^^%-^u;(m), and there- 



/..<hminfi^i4^. (4.19) 

m^oo |in|pi. 

From the bounded sequence {p{ql))q^n choose a convergent subsequence ip{qk^))keN 
and set p* := hmfc_^ooP('?fcl) G [0) !]• Then conditions (4.8) hold for the sequences 
Qkl and s((7fcl) with p*, and therefore (4.9) yields 

loga;(gfcl) , . 

By the definition of liminf we have liminfm^oo ^"^[^"^^ — hmfc_>oo log^^^fei)^ Hence 
by (4.19) and (4.20) we obtain ha < Xd{p*)- ' □ 



Proposition 4.2 Let d G N. Then for each m G N'^ 

log#Ty(m) ^ ^ log #1^0 (m) 
|m|pr - - 

log#l^(m) ^-^^ log#T^o(m) 2 
|m|pr ~ ~ |m|pr 

These upper and lower bounds converge to and hd, respectively, hence the latter 
are computable. 

Proof. The upper bounds follow from the general theory of NNSOFT (1.2), and 
their convergence from (1.1). For the lower bounds, let A: G N and consider the box 
(km). It can be decomposed into k'^ shifted copies of (m). Hence 

#Wo{kin) > #Wo{nif, #Wo{kni) > #Wo{nif . 

Sending A; to oo and using (4.6) and (4.7), we deduce the lower bounds as well as 
their convergence. □ 

We conclude this section by computing the various quantities in question for 
d = 1 and illustrating Theorem 4.1 for that case, where everything can be found 
explicitly. ^Wo(m) is the number of monomer-dimer tilings of (m). Clearly it 
satisfies #Wo(l) = 1, #Wo{2) = 2 and #Wo{m) = #Wo{m - 1) + #Woim - 2) for 

m > 3. It follows that #Wo{m) = F^+i, where Fm = ^ (^^)"* - 7f {^^Y 
are the Fibonacci numbers. #Wpcr(?7T') is the number of monomer-dimer tilings 
of T{m), and it satisfies #Wpor(l) = !> #^per(2) = 3 (one monomer tiling and 
two dimer tilings), and #Wpcr('7i) = i^WQ{m) + #Wo(m — 2) for m > 3 (the 
second term counting the tilings with a dimer occupying 1 and m). It follows that 

#VFper(m) = Fm+i + = ^m,, where Lm = [^^Y + (^^)"' are the Lucas 
numbers. if^W{m) is the number of monomer-dimer covers of (m), where a dimer 
may protrude from 1 to , or from m to m + 1. It satisfies ^W{1) = 3, ^W{2) = 5 



14 



and #VF(m) = #Wo(m) + 2#Wo(?w - 1) + #Wo("i - 2) for m > 3 (the three 
terms representing covers with zero, one, or two protruding dimers, respectively). 

It follows that #T^(m) = + 2F^ = + ^) (^)" +(1-75) ' 
From these values we see that ^°s#W{m) iog#VKpcr(m) iog#Wo(m) (,Qj^ygj.p.g 

m ' m m ^ 

hi = log ^-'y^, in accordance with (4.6). 

To determine Ai(p), it is enough to consider rational p £ [0, 1] by continuity, and 
then only n £ N such that s = ^ G N, and send such n to 00, by(4.8)-(4.9). Then 
7^Wo(s,n) is the number of linear arrangements of s dimers and n — 2s monomers, 
which is equal to ("7'*)- application of Stirling's approximation then gives [16] 

Xi{p) = lim ilog (^^ 2) 

n^oo n \ ^ 

1 - I) log (1 - I) - I log I - (1 - p) log(l - p). 

We see that Ai(0) = and Ai(l) = = /ii in accordance with (4.10) and (4.11). It 
is straightforward to verify that 

max Xi{p) = Ai (1 - = log ^ \^ = hi, 

in accordance with (4.12). 



5 Lower Bounds for Monomer-Dimer Entropy with Dimer 
Density p 

For an m X n matrix A, denote by perm^ A the sum of the permanents of all s x s 
submatrices of A. For a graph G, a matching is a set of vertex-disjoint edges, and 
W{G, s) denotes the set of all matchings of size s in G, which can be regarded as 
covers of the vertex set V{G) of G by s dimers (edges) and |V^(G)| — 2s monomers 
(vertices). If G is a bipartite graph with color classes (m) and (n), its incidence 
matrix is the mxn 0-1 matrix A = A{G) such that aij = 1 if and only if {i,j} is an 
edge of G. In that case it is immediate that ^W{G,s) = perm^ 74(G). A bipartite 
graph G is said to be r-regular if each vertex of G has degree r, equivalently A{G) 
has all row sums and column sums equal to r, so that ^A{G) is doubly- stochastic 
(a nonnegative matrix with all row sums and column sums equal to 1, necessarily a 
square matrix). 

Theorem 5.1 Let G be an r-regular bipartite graph with n vertices in each color 
class. Then 

WG,,)>(»)^!(r)'. (5.1) 

Proof. A result of the first author [11] states that if S is a doubly-stochastic 
nxn matrix, then perm^ B > perm^ J„, where J„ is the nxn matrix with all entries 
equal to ^. Since ^A{G) is doubly-stochastic, perm^ ^^(G) = ^pernig A{G) and 
perm^ J„ = (g)^^, the result follows. □ 
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The recent result of Schrijver [31] improves this lower bound for the case s = n 
if r is constant and n tends to infinity: under the assumptions of Theorem 5.1 

#WiG,n)>( ^'~}J~' ) . (5.2) 



It would be of interest to similarly improve the lower bound of Theorem 5.1 in the 
interesting range n large and s/n > r > (see below). 

In a recent paper [32], Wanless gives an alternative lower bound to (5.1), namely 
^W{G,s) > (") ^ — ^ • It turns out that except for ^ close to 1, the bound 
(5.1) is better. 

Theorem 5.2 Let d & N, p £ [0, 1] and recall the definition of \d{p), the 
monomer- dimer entropy with dimer density p, given by (4.8)-(4-9). Then 

>^d{p) > ^(-plogp- 2(1 -p)log(l -p) +plog2d-p). (5.3) 

Furthermore, the dimer entropy and monomer- dimer entropy h^ satisfy 

hd = Ad(l) > ^((2d - 1) log(2(i - 1) - {2d - 2) log 2d), (5.4) 

hd > ^i-pid) logp{d) - 2(1 - p{d)) log(l - p{d)) +p{d) log 2d - p{d)), (5.5) 
where 

Ad + l- yU+l 

p{d) = . (5.6) 

Proof. Let m = (mi, . . . , nid) G N'^ and assume that mi, . . . , m^ are all even. 
Let G be the adjacency graph of T(m). That is, the color classes of G are the sets 
{i G r(m) : zi + • • • + irf even} and {j G r(m) : ji + • • • + id odd}, and {i, j} is an 
edge of G if and only if i and j are neighbors on T(m), i.e., j = iitefc for some k € {d), 
where the addition is the standard addition in the group (Z/miZ) x • • • x (Z/m^Z). 
Then G is a 2d-regular bipartite graph on 2n = |m|pr vertices, and W{G,s) is the 
set Wper(in, s) of monomer-dimer covers of T(m) having exactly s dimers. Theorem 
5.1 yields that #Wper(m, s) > (")^s! {^Y ■ There is an injection / from Wper(m, s) 
to Wo(m + l,s), the set of monomer-dimer tilings of (m + 1) having exactly s 
dimers. If c G VFpcr(m, s), then /(c) is obtained from c by replacing each dimer 
in c occupying the points i = {ii,. . . ,id) and j = i + e^ such that ik = ruk and 
jk = by a dimer occupying the points i and (ii, . . . ,ik-i,mk + l,ik+i, ■ ■ ■ ,id)- 
Therefore #Wo(m + l,s) > (^)%! (f Let (mg)geN C N'^ and {sq)gm C N be 
sequences such that all the coordinates of each rrig = oo and 

limq_^oo \^rY~ ~ ^1 ~ ^"^2'^" ■ Then conditions (4.8) hold, and consequently 

(4.9) does. Therefore 

log #14^0 (mg,Sg) _ log#Wo(m5 + l,Sg) 



Xd{p) = lim j = lim 



> lim ^ "^^ — — — / j (piiy, 

q^oa 2nq n^oo 2n \pn J \ n 



pn 
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Manipulating the limit in the right-hand side of the inequality above and using the 
equality lim^^oo f (logr! — logr^) = —1, we deduce the inequality (5.3). 

Let (mg)ggN again satisfy the assumptions that all the coordinates of each lUg 
are even and limg_»oo = co, but this time set Sg = ng = li^^aifii. Using the 
inequality (5.2) for #Wpcr{mq,ng) and (4.11), we deduce the inequality (5.4). 

To prove (5.5), we use (4.12). We easily verify that the right-hand side of (5.3) is 
a strictly concave function of p in [0, 1], and p{d) given in (5.6) is its unique critical 
point in that interval, hence its maximizing point there. □ 

For d = 2,3, inequality (5.5) yields 

h2 > 0.6358077435 (5.7) 
/i3 > 0.7652789557. (5.8) 

For d = 3, inequality (5.4) yields /13 > 0.440075842, which is the best known lower 
bound. 



6 Upper and Lower Bounds on hd and hd Using Spectral 
Radii 

For d S N, X C (d) and m E N*^, we denote by (m/<) the projection of (m) on the 
coordinates with indices in K. Let VFper,_ft'(m), respectively VFper,A'(m), be the set 
of monomer-dimer covers, respectively dimer covers, of T{mK) x (m^^^\/^). Thus 
Wper^(d)(m) = Wper(m) and Wpor,(d)(m) = VFper(ni). Note that by the isotropy 

of our r, #Wper,E'(ni) and #Wper,_ft'(m) are invariant under permutations of the 
components of m if undergoes a corresponding change. 

In order to analyze Wpgj. |^|(m), we focus on the dimers in the cover lying along 
the direction e^. More precisely, with m' = (nii, . . . ,md-i), we consider (m') x 
Tirrifi) as consisting of irid levels isomorphic to (m'). A subset S of the points in 
level q is covered by dimers joining levels q — 1 and q (with level understood as 
level md); a subset T disjoint from S is covered by dimers joining levels q and q + 1 
(with level md + 1 understood as level 1); and the remainder U of level q is covered 
by monomers and dimers lying entirely within level q. We are interested in counting 
the coverings of U subject to various restrictions. With that in mind, for m' G N*^^^ 
we define an undirected graph G(m') whose vertices are the subsets of (m') in which 
subsets S and T are adjacent if and only if 5 n T = 0. When S n T = we also 
define, using U = (m') \ (5 U T), 

asT = number of monomer-dimer tilings of U 

bsT = number of monomer-dimer tilings of U viewed as a subset of T(m') 
PST = number of monomer-dimer covers of U, viewed as a subset of 
r(mi) X ((777-2, • • • jTrid-i)), each monomer within U, and each 

dimer meeting U but not S L)T. 
csT = number of monomer-dimer covers of U, each monomer within U, 

and each dimer meeting U but not S UT. 
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Thus in the tihngs/covers counted by asr, bsT, PST, cst, each monomer hes within 
U and each dimer meets U but not S UT. In asT, each dimer occupies two points 
of U that are adjacent in (m'). In bsr, each dimer occupies two points of U that 
are adjacent in T(m'), so is allowed to "wrap around". In psT, the dimers in the 
direction of ei are allowed to "wrap around" and the other dimers are allowed to 
"protrude out" of ((m2, • • • , In cst, tbe dimers may "protrude" out of 

(m'). Therefore asr < ^5T < PST < cst- By definition, iiU = 0, then asr = ^5T = 
PST = Cst = 1- Notice that when d = 2, there is no distinction between bsT and 

PST- 

We define the matrices A(m'), B{m'), i-*(m'), C(m') with rows and columns 
indexed by subsets of (m') as follows: 



A{in')sT 
B{m')sT 
P{m')sT 
C{m')sT 

Thus A{m'), B{m'), P(m'), C(m') are symmetric matrices — here is the "hidden 
symmetry" referred to in Section 4 — of integers satisfying < A{in') < B{m') < 
P(m') < C(m') (where the inequalities indicate componentwise comparisons). We 
denote by a(m'), /3(m'), 7r(m'), 7(m') their spectral radii, respectively, and conse- 
quently a(m') < P{m') < 7r(m') < 7(m'). 

In an analogous way, we define asT, bsT, PST, csT, where there are no monomers 
in the tilings and covers, the matrices A{m'), B{m'), P(m'), C(m') and their spec- 
tral radii 5(m'), /3(m'), 7r(m'), 7(m'). 

Each of these eight symmetric matrices can be considered as the adjacency ma- 
trix of an undirected multigraph, where the multiplicity of an edge is the correspond- 
ing matrix entry. This multigraph is a weighted version of G(m'). If the multigraph 
is bipartite, we say that the matrix is bipartite; if the multigraph is connected, we 
say that the matrix is irreducible; if the multigraph is disconnected, we say that 
the matrix is a direct sum; if the multigraph is connected and the greatest common 
divisor of the lengths of all its closed walks is 1, we say that the matrix is primi- 
tive, equivalently for sufficiently high powers of the matrix, all entries are strictly 
positive. 

Proposition 6.1 iei 2 < d G N and m = (m',md) G N'^. Then 

(a) tr A(m')™''* is the number of monomer-dimer tilings of (m') x T{md) and 
trA(m')"^'* is the number of dimer tilings of {m') x T{m^; 

(b) trB(m')™'* = #VFpcr(m) and tr 5(m')"'' = #Wpcr(m); 



asT 


if 5 n T = 





if 5 n T / 


bsT 


if 5 n T = 





if 5 n T / 


PST 


if s n T = 





if s n r / 


CST 


if 5 n T = 





if 5 n r / 
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(c) tTPiuiT" = #W^pcr,{i,d}(m) and trP(m')'"'' = #T^per,{i,d}(m); 

(d) tiCiuiT" = #Wper,{d}{m) andtiCimT' = #^^pcr,{d}(m); 

(e) for rud > 2, if column vector x = {xs)sc{in') is given by xs = bgijf, then 
x^S(m')"'''"^x = #Wpci.,(d_i)(m), if vector y is given by ys = Cs0, then 
y C(m')™'*~^y = #VF(m), and if z = {zs)sc{in') is given by zs = Psd, then 
z P(m')™'*^^z = #Wpf,j.,{i}(m); if column vector 3c = (Ss)sc(m') ^■^ sx^'en by 
Xs = bs9, then x^5(m')"d-25 ^ #VFpcr,(d_i) (m), i/y is given by ys = Cg^, 
then y C(m')™'*~^y = #VF(m), and if vector z is given by zs = Ps^, then 
i^P(m')"^'*-2i=#W?per,{i}(m); 

(f) the matrices ^(m'), B{m!), P(m'), C(m') are primitive; 

(g) if |m'|pr is odd, then A(m'), i3(m') are bipartite, otherwise they are direct 
sums. 

Proof. We begin with proving the first part of (b), its second part and (a), 
(c), (d) and (e) being similar. Assume first that = 1, and let (j) ^ W^por(m)- 
Since (j) can be extended periodically in the direction of with period 1, it can be 
viewed as an element of Wper(m'). Therefore 7^Wper(m) = #VFper(m'). We have 
trB(m') = X^5c(m') ^ss- Only the term 5 = contributes to the sum, and for this 
term we have U = (m') and 600 = #Wper(m'). Hence tri?(m') = #VFper(ni'). Now 
assume that rrid > 1, and consider any closed path 5i, 52, . . . , 5^^, Si of length 
in G(m'). For each p' G Sq place a dimer occupying the points (p', q) and (p', q + l) 
(with irid+l wrapping around to 1). We want to extend these dimers to a monomer- 
dimer tiling of T(m') x T(mrf) = T(m), i.e., to a member of Wper(m), by monomers 
and by dimers not in the direction of e^, i.e., lying within the levels 1, . . . , rrid- The 
number of choices of such monomers and dimers to fill the remainder of level q 
is given by 6s,_iSq) and so the number of extensions to a member of M^pcr(in) is 
bsiSibSiSa ■ • • bs^^_is„,^bs^^si- Conversely, each member of VFpcr(m) is obtained in 
this way. Hence #Wper(m) is the sum of all the products of the above form, namely 
trB(m')™''. 

To prove (f), we note that A{m') is irreducible, since whenever 5nT = 0, U can 
be tiled by monomers and therefore each subset of (m') is adjacent to in the graph 
of ^(m'). Furthermore, A{m') is primitive since the graph has a cycle of length 1 
from to 0. Since A{m') < B{m') < P(m') < C(m'), -B(m'), P(m') and C(m') 
are also primitive. 

To prove (g), let £, O denote the subsets of (m') with even and odd cardinality, 
respectively. If bsT > 0, then U can be tiled by dimers and so #U must be even. 
Therefore if |m'|pi. is odd, members of £ are adjacent only to members of O in the 
graph of -B(m'), and so that graph is bipartite; if |m'|pr is even, then members of £ 
are adjacent only to themselves, and the graph is disconnected. The same conclu- 
sions hold for A(m') since A(m') < B(m'). □ 
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Lemma 6.2 Let 2 < d e N and m' e N'^"^ Then 

(mO (6.1) 

= log/?(m') (6.2) 

log7r(m') (6.3) 

m') (6.4) 

(m') (6.5) 

= log ^(m') (6.6) 

log^(m') (6.7) 

m'). (6.8) 

Proof. From Part (a) of Proposition 6.1 we obtain ^W(){m',md) < tr^(m')™'*, 
and therefore 

log#H^(m>id) . 1- logtrl(m0'"'^ , 
limsup < Imisup = loga(m j. (6.9) 

m^j— >oo TTld m^j^oo "Id 

The equality in (6.9) follows from a characterization of p{M) for a square ma- 
trix M > 0, namely p{M) = limsup„^oo(tr M")~ (see for example Proposition 
10.3 of [13]). Since — log ^Wo{m' ,md) is subadditive in md, the first limsup in 
(6.9) can be replaced by a lim, which proves (6.5). Similar considerations prove 
limm^^oo '"'''^ < loga(m'). In order to prove the reverse inequality and 

thus (6.1), observe that each monomer-dimer tiling of (m') x T{md) extends to 
a monomer-dimer tiling in VFo(m'i"T'd + 1) (replace each dimer occupying (m',1) 
and (m',md) by a monomer occupying (m',1) and a dimer occupying (m',md) 
and (m',md + 1), and tile the rest with monomers). Hence #T4^(m',md + 1) > 
trA(m')™''' by Part (a) of Proposition 6.1. Therefore, since — log T^VFo(m', md) is 
subadditive in md and thus the limits below exist, we obtain 

j.^^ log#VFo(m',md) _ ^.^ log#VFo(m',md + l) 



md^oo md m^-»oo md 

logtr^(mO"'^ 
> limsup = loga(m j. 

m^j— >oo "id 

To prove (6.2), (6.3) and (6.4), we use another characterization of the spectral 
radius. A vector norm is a mapping || • || : Mn{C) — > M-|- taking complex matrices 
of order n to nonnegative reals such that ||M|| = only if M = 0, = |2;|||M|| 

for all z E C, and \\M + N\\ < \\M\\ + \\N\\. If aj > for all i,j G (n), then 
||M|| = Ylij is a vector norm. Proposition 10.1 of [13] states that if || • || is 
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a vector norm, then p{M) = limfc^oo ll^'^^lh- In particular, if M > and v is a 
column vector with positive entries, then p{M) = limfc_>oo(v^M'^v) . Applying this 
to M = B{m'), P(m'), C(m') and using Part (e) of Proposition 6.1 with v = x, z, y 
defined there proves (6.2), 6.3), (6.4). 

The proof of (6.6) is a little more complicated because the vector x in Part 
(e) of Proposition 6.1 is not strictly positive. Therefore we introduce the vector 
w with entries ws = max(l,X5). Then, by Part (e) of Proposition 6.1, we have 
#W^pcr,(d-i)(m) = x^i?(m')™''~^x < w^i3(m')"*''~^w. Therefore we obtain 

log#VFper,(d_i)(m) logw^5(m')'"''-2w , ~ . 
hm ^ < hm = logp(m ) 

(the first Hm above exists since log ^Wpor.(d-i) (m) is subadditive in md). On the 
other hand #H^per,(d-i) (m) > #l^per(m) = tr i?(m')™'* by Part (b) of Proposition 
6.1. Therefore 

,. log#M?per,(d-i)(m) logtr5(m')™'* . ~, 
iim > hmsup = iogp(m ). 

This proves (6.6). To prove (6.8), we show analogously that 

log#W>i) 
nm < log7(m j, 

and on the other hand, by Part (d) of Proposition 6.1, 

log#Ty(m) log#W^per,{d}(m) 
lim > hm sup ■ 



= hmsup = iog7(m j. 

The proof of (6.7) is similar. □ 



Proposition 6.3 Xe< 2 < d G N and m' G . Then 



Im'l 



pr I "J- \pr 



log7^^~^^loga^_ (6.11) 



\pr I"'- \pT 



Proof. The upper bounds follow from the general upper bounds in Proposition 
4.2 along with (6.4), (6.8). The lower bound in (6.10) follows similarly from the 
general lower bound in Proposition 4.2 along with (6.1). However, since (6.5) only 
gives a lower bound for loga(m'), we use a separate argument for the lower bound 
in (6.11) as follows. For g G N, ^W{qm! ^rud) is not smaller than the the number 
of dimer tilings of {qm') x T(md), which in turn is not smaller than the number 
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of dimer tilings of (m') x T{md) raised to the q'^ ^ power. Hence by Part (a) of 
Proposition 6.1 we have 



and so 



Therefore 



log i^W {qm' , md) ^ logtvA^'^ 
\{qin',md)\pr ~ |m'|pi.md ' 



r 1- log#VK((7m', m^) 1 logtrA™'* loga(m') 
hd= hm — >- — — hmsup = — — . 

q,ma^oo \[qin',md)\pr |m'|pr 5,md^oo rUd \pr 

□ 

Now we introduce the following notation. For m G N'^ and fc E (d), m^^ := 
(mi, . . . , mfc_i, . . . , rud) G N'^"^. As special cases we have the previous nota- 

tion m' = m~'^ and m~ = m~^. 

Proposition 6.4 Let m G N"', and assume that rud is even. Then each k E 
{d — 1) satisfies 

log/3(m~'^) log2 log/3(m~'=) 

\ — \ ^ \ — ZFt (6.12) 

|m|pi. mfc \m^^\pr 

log^(m~'^) ^ log2 ^ log/3(m~^') 

m 7 I nrt'^f^ I 



|m|pr ruk |m' 



Proof. We have 



/3(m~'^)'"'* < trB(m~'^)'"'* = #Wpcr(m) = tr S(m~*^)"''= 



' per 



< 21""^ lp■^/?(m~'=)"'^ 

where the first inequality follows since /?(m~'^) is one of the eigenvalues of i?(m~'^), 
which are all real, and nid is even, the next equality from Part (b) of Proposition 
6.1, the next equality from the same and the fact that ^VFpor(m) is invariant under 
coordinate permutations in m, and the last inequality from the fact that i?(m~'^) 
has 21"^ eigenvalues, all real, whose absolute values are at most /3(m~'^). Taking 
logarithms and dividing by |m|pr, we deduce (6.12). The inequality (6.13) is obtained 
in a similar way. □ 
We define 

— , ^ log #TVncr •rii("ii, m^) _ ^ ^ 



hd^iimi) := hm ^ N; /irf_i(0) := log2. 



Notice that for mi G N, hd-i{mi) is the same as /i(mi,r) defined in (3.1) when 
r is the d-digraph encoding the monomer-dimer covers. For this reason the limit 
hd-i{mi) exists, and similarly for hd-i{Tni). The following theorem is an analog of 
Theorem 3.1 and Theorem 3.4. 
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Theorem 6.5 Let 2 < d e 'N, p,r e 'N, q£'L+. Then 

hd~ii2r) ^^^^ hd^i{p+2q)-hd-i{2q) 

2r ~ ~ p 

hd-i{2r) ^-^^ hd-ijp + 2q) - hd^i{2q) 

2r ~ ~ p 

Let m' = (mi, . . . , md-i) G N"^"^ and assume that mi, . . . , m^-i are even. Then 

hd < (6.16) 

Proof. We have 

r log#Wo(m',"Zd) log5(m') log7(m') 
hd = lim j — < iim mi — j — < hm sup — j — 

m'.m^^oo m prm^ m'^oo m pr m'^oo m pr 

log#VF(m',md) ~ 
= hmsup ■ — ^1 = hd, 

where the first equahty follows from (4.7), the next inequality from (6.5), the next 
one from 5(m') < 7(m'), the next equality from (6.8), and the last equality again 
from (4.7). From this and 5(m') < /3(m') < 7(m') we obtain 

hd= hm l^g^^ j.^ Iog7(m0_ (g_,8) 

m'— »oo |ni Ipf m'— »oo |ni |pi- m'-^oo |in |pi- 

Similarly (and more simply) 

hd= hm ]^f^= li^ ) = hm (6.19) 

m'— >oo |m Ipr m'^oo |m |pi- m'— >oo |m [pi- 

First we prove (6.16). Let m' = (mi, . . . ,md-i) € N, mi, . . . ,md-i even, and 
let p = (pi, . . . ,Pd~i) G be arbitrary. Set 



n ■ ■ ■ 5 



mi = (pi, . . . ,pd^i,mi), m2 = {p2, ■ ■ ■ ,Pd,mi,m2), 

m^-i = {Pd,mi, . . . ,mrf_i). 

Then, using (6.12) with k = 1 d — 1 times, we obtain 

log/?(p) ^ log 2 ^ log/3(mj;) ^ log 2 ^ log 2 ^ log/3(m2) ^ 

\p\pr ~ Pi \l^l\pr ~ Pi P2 11^2 \pr ~ 

^ glog2 ^ log/3(mO 

Letting p ^ oo and using (6.19) for the left-hand side, we deduce (6.16). Similar 
arguments apply to deduce (6.17). 
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We now demonstrate the lower bound in (6.14). Let m~ G N'^"^, p G N, g G 
Assume first that g G N. Since 7(01^) = p(C(m~)) and C(m^) is symmetric, it 
follows as in the arguments for (3.8) that 

> tj^Cirnr)^ ^ #^,.,,,,(p + 2,,m-) 
^ - trC(m-)2'? #T^p,,,^i}(2cz,m-) 

Taking logarithms, dividing by | m \pr^ letting m — > cx), and using (6.19) and the 
definition of /i(i-i("^i)) we deduce the lower bound in (6.14) for the case g G N. If 
g = 0, we have to replace the denominators in (6.20) by tri = 2l™ lp'', and the 
lower bound in (6.14) is verified because hd-i{0) was defined to be log 2. The lower 
bound in (6.15) is proved similarly. 

We now prove the upper bound of (6.14). For each m' G N^~^ we have 

7(m')''" < trC(m')'^ = #VFper,{4(m', 2r) = #VFpe,,|i|(2r, m'), 

where the inequality above is true because the eigenvalues of the symmetric ma- 
trix C(m') are real and 7(m') is one of them, the first equality follows from Part 
(d) of Proposition 6.1, and the last equality from the invariance under coordinate 
permutations. Therefore 

Iog7(m0 ^ log#^per,{i}(2r,m0 
|m I 2t" I m I 

and letting m' — > 00, we deduce the upper bound of (6.14) by (6.19) and the defi- 
nition of hd~i{m\). Similarly we deduce the upper bound of (6.15). □ 

The following theorem supplies practical upper and lower bounds on 2- and 
3-dimensional monomer-dimer and dimer entropies. 

Theorem 6.6 Let p, r,t,u,v G N and q,s £ Then 

log/?(2r) ^ ^ log/?(p + 2g)-log/3(2g) ^ ^ 

2r ~ ~ p 

log/3(2r) ^ ~ ^ log^(p + 2g)-log^(2g) ~ 

^ ^ ^2 > , /^(O) = 2 

log /3(2r, 2t) log f3{p + 2q,u + 2s) - log p{p + 2q, 2s) log P{2q, 2v) 

^ /13 ^ 

Art up 2vp 

log /3(2r, 2t) ^~ ^ log (3{p + 2q,u + 2s) - log P{p + 2q, 2s) log j3{2q, 2v) 

> /13 > 

Art up 2vp 

(3{n, 0) = /3(0, n) = (3{n, 0) = ^(0, n) =2", n G N. 

Proof. The upper bounds in the above inequalities are the inequalities (6.16) 
and (6.17). We now show the lower bounds. Equations (6.2) and (6.6) for d = 2 
yield 

hi{mi) = log (3{mi), /ii(mi) = log/3(mi), mi G N. (6.21) 

Hence the lower bounds on /i2,/i2 follow immediately from the lower bounds in 
(6.14), (6.15), equation(6.21) and the equalities 7ii(0) = ^i(O) = log 2. 
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In order to establish the lower bounds on on h^jh^, we first establish lower and 
upper bounds on /i2("T'i) and h2{mi) in terms of /?(•, •) and •). The definition of 
/i2(mi) and /i2("ii) and equations (6.3) and (6.7) for d = 3 yield 

T t \ r logvr(m') r . x y log7f(m') 

h'2{n^i) = hm , ri2[mi) = lim , mi G (6.22) 

where m' = (mi, m2). Since P(m') is a nonnegative symmetric matrix with spectral 
radius 7r(m'), it follows as in (3.8) and using Part (c) of Proposition 6.1 that 

, trP(mO"+2^ _ #VFpe,,{l,3}(m^n + 2s) 
""^"^^ - trP{m'r^ #t^p,,,|i,3}(m',2s) ■ 

Here ti G N and s G Z+. When s = 0, trP(m')^'' = 2l'"'lp'', and so this is the value 
we use for ^Wpg^ ,[i 3}(m', 0). Take logarithms of this inequality, divide by m2 and 
send m2 to oo. Using (6.22) and (6.2) for = 3, we deduce that 

7^2(rr^i)>^°"^^"^^'^ + '^^"^°^^^^^''^\ miGN, (6.23) 

u 

where /?(mi,0) := 2^^. Similarly 

/^2(mi)> ^°g^^"^^'" + '^^"^°g^^"^^''^\ miGN, (6.24) 

u 

where /3(mi,0) := 2™i. For -y G N we have the inequality 7r(m')^*' < trP(m')^^ = 
#^per,{i,3}(i^') 2"^). Take logarithms of this inequality, divide by 2vm2 and send 
m2 to CO. Using (6.22) and (6.2) for d = 3, we deduce that for mi G N 

7^2(mi) < (6.25) 
2v 

Inequality (6.25) also holds for mi = since by definition h(0) = log 2 and /3(0, 2v) = 
2^". Similarly, for mi G Z+ 

him,) < (6.26) 
2v 

Now we can substitute the bounds (6.23) and (6.25) in the lower bound of (6.14) 
as appropriate from the signs in the numerator, and obtain the lower bound on /13 
as stated in the theorem, and similarly for ^3. □ 



7 Using Automorphism Subgroups to Reduce Compu- 
tations 

The matrix B(m') has order 2", where n = |m'|pi., and so has 4" entries. Since its 
{S, T) entries are positive precisely when S* n T = 0, its number of positive entries is 
^ (^)2"~* = 3". Hence it is sparse. However, already for m' = (4, 4) it has 4.3 • 10^ 
nonzero entries, and the computation of its spectral radius is infeasible for standard 
PC. Nevertheless, this computation can be reduced to computing the spectral radii 
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of a suitable nonnegative matrix whose order is the number of orbits of the action 
of an automorphism subgroup of B{m'). This usage of automorphisms is also used 
in [7] and [26]. 

Recall that given an NxN complex-valued matrix A = {aij)i , its automorphism 
group is the subgroup of the symmetric group Sn on {N) defined by 

Aut(A) := {it £ Sn ■ a7r(j)7r(i) = o-ij for all i,j £ {N)}. (7.1) 

Let ^ be a subgroup of Aut(^). The action of Q partitions {N) into minimal 
invariant subsets called orbits. We denote by O := {N)/Q the orbit space (set of 
orbits), and by Greek letters a,(3, . . . its members. We have 

'^o.ij = '^(^TT{i)TT{j) ='^a^{i)k, a,l3£0, i £ a, tt £ Sn, (7.2) 

which means that for given a, (3 £ O, the sum Tij^paij is the same for all i £ a. Let 
M = and define the M x M matrix A = {aafi)a,i3&o by 

aa/3 = ^ajj, i£a. (7.3) 

ie/3 

This is a valid definition by (7.2). The following proposition is known, and we prove 
it for completeness. 

Proposition 7.1 Let A = {aij)i he a complex-valued matrix. Let Q he a suh- 
group of Aut A, O its orhit space, and M = i^O. Let A be the induced M x M 
complex-valued matrix given hy (7.3). Then the spectrum (set of eigenvalues) of 
A, spec(A), is a subset o/spec(^), and in particular p{A) < p{A). If A is a real- 
valued nonnegative matrix, then p{A) = p{A). If A is real and symmetric, then A 
is symmetric with respect to an appropriate inner product on M*^, and in particular 
spec(yl) is real and A is diagonalizable. 

Proof. Let n^r be the group of N x N permutation matrices. Let l : Sn 
be the standard representation of Sn- That is '-(vr)(xi)jg^Ar^ = (x7r(i))ig(Ar). Let 

A" := {x G : t(7r)(x) = x for all tt £ G} 

= {{xi)i(z(^N) G : = Xi for ah i £ {N),Tr £ Q] 

be the subspace of vectors that are constant on each orbit of Q. Then X C is 
the largest subspace of on which i{Q) acts trivially (as the identity operator). 
Clearly, X is isomorphic to C*'^. Indeed, each -x. = {xi) £ X induces a unique vector 
x := {xa)aeo S C^, where Xq, = Xi for any i £ a. Conversely, each y £ induces 
a unique £ X such that y = x. Next, we observe that X is an invariant subspace 
of A. Indeed, for each x = (xj) £ X and vr G ^ we have for all i £ (N) 

N N N 

(ylx)j = C'ij^j — 5^ ^7r(i)7r(j)^7r(j) = ^Trfpfc^fc = {^'^)iT{i), 

j = l j = l k = l 

which means that Ax £ X . Moreover, ii 'x. £ X and x = (xq) £ C^'^ is defined 
as above, then for any i G a we have (Ax)j = J2f3eo^cef3Xf3, and consequently 
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Ax = j4x. This means that the action of j4|;^:' is isomorphic the the action of A on 
C*^. In particular, 

spec(y4) = spec^Alx) ^ spec(74), 

and therefore 

p{A)<piA). 

Assume now that A is nonnegative. Then by the Perron-Frobenius theorem p{A) S 
spec(A), and A has an eigenvector x belonging to p{A). Since each vr G Aut(A) 
satisfies Al{tt) = l{tt)A, it follows that l{tt)x is also an eigenvector of A belonging 
to p{A). Hence X^7reAut{A) '-i'^)^ G ^ is an eigenvector of A belonging to p{A). 
Therefore p{A) S spec(A|;v:') = spec(yl). It follows that p{A) = p{A). 

Finally assume that A is a real symmetric matrix. That is {Ax,y) = (x, Ay), 
where (x, y) = y^x is the standard inner product in M^. For each a G O, let Wa 
be the cardinality of the orbit a. In we define the inner product 

(x,y) := '^WaXaVa- (7.4) 

Then all x, y G A" satisfy (x, y) = (x, y). Hence (Ax, y) = (x. Ay), i.e., A is sym- 
metric (self adjoint) with respect to the inner product (•, •) in M*^. In particular, A 
has real eigenvalues and is similar to a diagonal matrix. □ 

We shall now briefly mention the power method for computing p{A) where A is 
a nonnegative symmetric matrix of order N, and a variant of it that works on A of 
order M, which we used in our computations. 

Proposition 7.2 Let A be a nonnegative symmetric matrix of order N . Choose 
a scalar r > and a positive vector xq = (xo,i, . . . , xq^n) ■ For each m G N, let 

Xm = iXm,l Xm,Nf ■= {A + rI)Xm-l 

Im ■■= mm 

* Xm—l,i 

Xm,i 

Um ■= max 

* Xm—l,i 
(XjT) J ^m— 1 ) 

\^m— 1 1 Xjrj— 1 j 

Then Im is nondecreasing and is nonincreasing in m, 

lm<rm< p{A) + r < Urn, m G N 

lim Im = lim Um = p{A) + r, 

m— >oo m— >oo 

and Xral \J (xrri) Xm) convcrges to an eigenvector of A belonging to p(A). 

Furthermore, with the notation of Proposition 7.1, if we choose the vector xq to 
he in X , i.e., ifxQ is constant on each orbit of Q, then each m G N the vector x^ is 
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also in X (so is defined), 



X 



-m 



{A + r7)X^_i 



m 



mm — 



m—l,a 



U, 



•m 



max — 



m 




and Xm/ 



X 



Xm) converges to an eigenvector of A belonging to p{A) = p{A). 



For m' G N'^"^, let GT(in') be the adjacency graph of the elements of the torus 
T(m'). The automorphisms of Gr(m') act as automorphisms of the symmetric non- 
negative matrices B(m') and B{in'). In view of Proposition 7.2, in order to compute 
the spectral radii /?(m') and /?(m'), it is advantageous to use large automorphism 
subgroups of GT(m'). The rigid motions of the box (m') and of the torus T(m') 
are automorphisms of Gr(m'). 

The rigid motions of (m') contain the reflections across the hyperplanes = 
"^''^^ ^ k £ (d—l), which commute with each other, and the allowable transpositions 
exchanging Xj and Xj in case = ruj. Thus if m' = ml, m >2, then the group of 
rigid motions of the cube (m') contains a subgroup of order 2'^~^((i— 1)!. For d = 2,3, 
which is our main focus in this paper, the reflections and allowable transpositions 
generate all the rigid motions of (m'). 

The rigid motions of T(m') contain, in addition to the rigid motions of (m'), the 
unit translations x x + e^, k £ (d—l). The unit translations generate the group 
of translations, an Abelian group isomorphic to (Z/miZ) x • • • x (Z/m^-iZ) of order 
|m'|pr- We call the group generated by the reflections, the allowable transpositions 
and the unit translations the group of rigid motions ofT{m'). Note that for T(2) 
the reflection coincides with the unit translation, and similarly for r(m'), if = 2 
then the reflection across = | coincides with the unit translation x i— > x + e^. 
We are aware of additional automorphisms of Gr(m') if at least two components 
of m' are equal to 4: observe that Gj'(4) is isomorphic to Gr(2,2), since both are 
4-cycles. Therefore Gt(4, 4) is isomorphic to Gt(2,2,2,2), and its automorphism 
group has order at least 2^ • 4! = 384, whereas the group of rigid motions of T(4, 4) 
has order 2^ • 2 • 4^ = 128. Similar results hold for d > 3. 

The following proposition is straightforward: 

Proposition 7.3 Let Pi, . . . , C (n) x (n) and P = (Pi, . . . , P^). Let m = 

(?7T,i, . . . , mrf) G N'^ and m' = (mi, . . . , ?7i^„i), and consider the transfer digraph 
0(i(m') between members o/Wpcr(m') with respect to Td- Then the group of trans- 
lations o/T(m') acts a subgroup of automorphisms o/0£;(m'). If for some k G {d—l) 
Pfc is symmetric, then the reflection across the hyperplane xj. = SUd— Q^ig 
automorphism of ©^(m'). If for some p,q £ (d—l) nip = nig and Tp = Tg, then 
the transposition exchanging Xp and Xg acts as an automorphism o/©rf(m'). 

Corollary 7.4 Let Pi, . . . ,Pfi C (n) x (n), P = (Pi, . . . ,Pd), and assume that 
Pi, . . . ,Pd_i are symmetric. Let m = (mi, . . . ,m^) G N'^ and m' = (mi, . . . ,md-i), 
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mi 




log/3(mi) 


log/3(-mi) 
m-] 


4 


6 


2.6532941163 


.66332352908 


5 


8 


3.3135066910 


.66270133821 


6 


13 


3.9769139475 


.66281899125 


7 


18 


4.6395628723 


.66279469604 


8 


30 


5.3023993987 


.66279992338 


9 


46 


5.9651887945 


.66279875494 


10 


78 


6.6279902386 


.66279902386 


11 


126 


7.2907885674 


.66279896067 


12 


224 


7.9535877093 


.66279897578 


13 


380 


8.6163866375 


.66279897212 


14 


687 


9.2791856222 


.66279897301 


15 


1224 


9.9419845918 


.66279897279 


16 


2250 


10.60478356551861 


.662798972844913 



17 4112 G (11.26758254, 11.26758315) G (.6627989729, .6627990088) 

Table 1: Spectral radii for /i2 

and assume that for all p,q £ {d — 1), Fp = Tq if rup = niq. Then the automor- 
phism subgroup of Gt{t^') described above acts as an automorphism subgroup of the 
transfer digraph 0d(m'). 

As an example, consider the upper and lower bounds given by (3.9). The pa- 
rameter 92{m) appearing there is the spectral radius of the matrix B{m) defined 
in Section 6, which has an automorphism subgroup of order 2m, isomorphic to the 
group of rigid motions of T{m), if m > 2. -6(15) is 2^^ x 2^^, but as we shall see, 
-6(15) is 1224 X 1224, which makes the computation of its spectral radius feasible 
on a regular desktop computer. 

These observations are our main keys in finding good upper and lower bounds for 
/i2 and /13. We point out that [7] was the first work that used these automorphisms 
of B{m') to help obtain a good upper bound for h^, which was later improved in 
[26] by similar methods. 

8 Numerical Results for Monomer-Dimer Entropy in 
Two and Three Dimensions 

Our results are based on Theorem 6.6, and we compute the spectral radii appearing 
there using Propositions 7.1 and 7.2, and the automorphism subgroups described in 
Section 7. We first consider the two-dimensional monomer-dimer entropy. Recall 
that j3{mi) is the spectral radius oiB{mi). Table 1 lists log/3(mi), ^^^^^p^, and the 
number ^0{mi) of orbits of the torus T(mi) under the action of the group of rigid 
motions of T(mi). The computation of log/3(17) was interrupted, and the table 
indicates the best interval in which we can locate it. We notice that the sequence 
^^^^2ll is decreasing for r = 2, . . . , 8. Hence /i2 < iHlMM = .662798972844913 is 
the best upper bound for /i2 from our data. The best lower bound for /i2 from our 
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log/3(mi) 
mi 
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39Q9 








980Q 




13 


1 8437Q7937 


3073 


7 


-LO 


^ ^ij^ 




Q 

KJ 


30 


9 400849903 


3001 


q 




9 'iQ4837310 


9883 


1 n 


7fi 
( 


z.yDyooyzo < 


./yoy 


11 


126 


3.183303939 


.2894 


12 


224 


3.543130579 


.2953 


13 


380 


3.770113562 


.2900 


14 


687 


4.119721251 


.2943 


15 


1224 


4.355934472 


.2904 



Table 2: Spectral radii for /i2 



data is /12 > > .66279897. This improves the lower bound (5.7) 

from permanents by more than 4%. Hence we obtain the value 

h2 = .66279897, (8.1) 

correct to 8 decimal digits. We also notice that the sequence ^"^2j'+i~^^ increasing 
for j = 2, . . . , 8. Suppose that this sequence were increasing for all values of j. Since 
linij^oo ^°^2j+^^^ = ^2 by (6.19), it would follow that /i2 > Mil) > .6627989729. 
The last digit of this bound is too high, as seen by comparison with our best upper 
bound, probably caused by roundoff errors in the interrupted computation, but 
enables us to state that the above hypothesis would gives the value ^2 = .6627989728 
correct to 10 digits, consistent with the one found by Baxter [1] (his value of /12 is 
accurate to 8 digits, as can be seen by evaluating log ^ for s = 1 in his Table II 
and varying the last digit of the tabulated ^). Since the lower bound (5.7) for /12 is 
quite close to the correct value of /i2, it is reasonable to assume that the value p*, 
for which X2ip*) = ^2, is fairly close to p(2) = ^^'^ ~ 0.6096118 (according to [1], 
p* = 0.63812311.). 

As a check. Table 2 gives /3(mi), the spectral radius of B{mi), yielding lower 
and upper bounds for the known entropy /12 = 0.29156090 .... Again, the sequence 
log ^{2r) dggpgg^ggg fQj, J. — 2, ... ,7 and the sequence ^"^2j+]^^^ increases for j = 

2, ... ,7. Thus the best upper bound on /12 from our data is i^i^M = .2943, which 
is larger by 0.9% than the true value. The best lower bound is ^°g/^(^^)~*°g/^(^2) _ 
0.2883, which is smaller by 1.1% than the true value. We notice that = 

.2905 < /i2, consistent with the assumed fact that ^"^2j+]^^^ increases for all j. 

We now consider the three-dimensional monomer-dimer entropy h^. Recall 
that /3(mi,m2) = P{m2,mi) is the spectral radius of B{mi,m2). Table 3 gives 
log /3(mi, 7722), ^"^mTm'r^^^ ' number #0(7771,7712) of orbits of the torus 

T{mi, 7712) under the action of the group of rigid motions of T(7ni, 7772). (In the case 
(7771,7^2) = (4,4), we recall that the group of rigid motion of T(2, 2) has order 128, 
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(mi,m2) 


#0(7711,7712) 


log /3(77ii, 7712) 


log/3(mi,m2) 


(2,2) 


6 


3.224405658 


0.8061014145 


(3,2) 


13 


4.768958913 


0.7948264855 


(4,2) 


34 


6.367778959 


0.7959723699 


(5,2) 


78 


7.958105292 


0.7958105292 


(6,2) 


237 


9.550024542 


0.7958353785 


(7,2) 


687 


11.14163679 


0.7958311993 


(8,2) 


2299 


12.73331093 


0.7958319331 


(3,3) 


25 


7.057039652 


0.7841155169 


(4,3) 


158 


9.421594940 


0.7851329117 


(5,3) 


708 


11.77517604 


0.7850117360 


(4,4) 


805 


12.57923752 


0.7862023450 



Table 3: Spectral radii for /i3 



(7711,7772) 


#0(7771,7772) 


log /3(r7ii, 7772) 


log /3(mi ,7712) 
mi 


(2,2) 


6 


2.292431670 


0.5731079175 


(3,2) 


13 


3.068671222 


0.5114452037 


(4,2) 


34 


4.151763891 


0.5189704864 


(5,2) 


78 


5.119835223 


0.5119835223 


(6,2) 


237 


6.161467494 


0.5134556245 


(7,2) 


687 


7.168058989 


0.5120042135 


(3,3) 


25 


3.938705096 


0.4376338996 


(4,3) 


158 


5.365527945 


0.4471273287 


(5,3) 


708 


6.635849120 


0.4423899413 


(4,4) 


805 


7.409698288 


0.4631061430 


(6,3) 


4236 


7.97716207 


0.443175671 


(6,4) 


184854 


10.98112634 


0.4575469308 




Table 4: Spectral radii for h 


3 



and it turns out to have 805 orbits. We also did the computations with the larger au- 
tomorphism subgroup of Gt(4, 4) of order 384 discussed in Section 7, which turns out 
to have 402 orbits. Both computations gave the same value of /3(4, 4).) Recall that 
^3 < '"^^irT'^*'* i hence the best upper bound for /13 from our data is ^^"^'^^ = 
0.7862023450. The best lower bound is '°g/3(3,5M°g/3(3,4) _ i°g^(2,8) ^ .761917234. 
It turns out that the permanent lower bound (5.8) is better: > .7652789557. 
Of course, had we computed f3{mi,m2) for larger ttii and 777-2, we would eventually 
improve the permanent lower bound. Thus, the best estimates we have are 

.7652789557 < < .7862023450. (8.2) 

Table 4 lists /3(777i, 7772), the spectral radius of -6(7771,7772), which give bounds for 
h^. The entry (7771,7772) = (6,4) is taken from [26], which took advantage of the 
fact that the matrix of order 184854 is a direct sum of 3 matrices. The best upper 

bound for is = 0.4575469308, which was reported in [26] . The best lower 



31 



bound from the data is given by '°g/3(4.6Mog/3(4,4) _ iog|^ ^ .3794013885, which is 
a weak lower bound. The best lower bound for is given by (5.4): h^, > 0.4400758. 

We now compare our results for /12 with the results of [19]. On page 342, 
Hammersley and Menon tabulate estimates of X2{p) computed by the Monte Carlo 
method in increments of 0.05 for < p < 1. The maximal value in their table is 
.6676 for p = 0.65. They state "There are reasons for believing that this Monte 
Carlo estimate has a small negative bias, probably 1% or 2% too low". However, 
since A2(p) < ^2 = .66279897, the Monte Carlo estimate for A2(0.65) is at least 0.7% 
higher than the true value. 

We conclude with a comparison of several lower bounds for the monomer-dimer 
entropy with dimer density p, Xd{p), for d = 2, 3. Hammersley and Mennon [19] give 
a lower bound for Xd{p), graphed and tabulated in increments of 0.05 for < p < 1. 
Bondy and Welsh [5] give another lower bound for \d{p), which depends on the dimer 
entropy Ad(l) and increases with it. Since A3(l) is known only through upper and 
lower bounds, the bound of [5] improves each time a better lower bound for Arf(l) 
is found. We computed the lower bound of [5] for Xsip) using the best available 
lower bound A3(l) = /13 > 0.4400758. Hammersley and Mennon too tabulated and 
graphed the bound of [5] for Xs{p), but at the time the available lower bound for 
A3(l) was weaker. Figures 2 and 3 illustrate the lower bounds for Xdip), d = 2,3, due 
to [19], [5], and Theorem 5.2. Figure 2 also illustrates the Monte Carlo estimates 
of [19]. It is seen that except for very high p, the best lower bound is given by 
Theorem 5.2. (As pointed out above, (8.1) implies that the Monte Carlo estimates 
above the line y = /i2 are over estimates). We also include in the figure estimates of 
X2{p) obtained from the heuristic computations of Baxter [1]. One can obtain from 
the lower bound of [32] a corresponding lower bound for Arf(p). It turns out that 
for d = 2,3, this bound is dominated by the maximum of the lower bound given by 
Theorem 5.2 and the lower bound of [5]. 
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h2 




Figure 2: Lower bounds and estimates for A2(p). HM is the lower bound of [19], 
BW is the lower bound of [5], FP is the lower bound of Theorem 5.2, MC is the 
Monte Carlo estimate of [19], B is the estimate from [1], and h2 is the true value of 
/i2 = max A2(p). 



0.8 T hSHigh 




P 

Figure 3: Lower bounds for A3(p). HM is the lower bound of [19], BW is the lower 
bound of [5], FP is the lower bound of Theorem 5.2, h3Low and h3High are the best 
lower and upper bounds for /13 = max A3 (p). 



33 



References 

[1] R.J. Baxter, Dimers on a rectangular lattice, J. Math. Phys. 9 (1968), 650- 
654. 

[2] R.J. Baxter, Hard hexagons: exact solution, J. Phys. A: Math. Gen. 13 
(1980), L61-L70. 

[3] I. Beiclil and F. Sullivan, Approximating the permanent via importance 
sampling with application to the dimer cover problem, J. Comp. Phys. 149 
(1999), 128-147. 

[4] R. Berger, The undecidability of the domino problem, Mem. Amer. Math. 
Soc. 66 (1966). 

[5] J. A. Bondy and D.J. A. Welsh, A note on the monomer dimer problem, 
Proc. Camb. Phil. Soc. 62 (1966), 503-505. 

[6] N.J. Calkin and H.S. Will, The number of independent sets in a grid graph, 
SIAM J. Discr. Math. 11 (1998), 54-60. 

[7] M. Ciucu, An improved upper bound for the 3-dimensional dimer problem, 
Duke Math. J. 94 (1998), 1-11. 

[8] M.E. Fisher, Statistical mechanics of dimers on a plane lattice, Phys. Rev. 
124 (1961), 1664-1672. 

[9] S. Forschhammer and J. Justesen, Entropy bounds for constrained two- 
dimensional random fileds, IEEE Trans. Info. Theory 45 (1999), 118-127. 

[10] R.H. Fowler and G.S. Rushbrooke, Statistical theory of perfect solutions. 
Trans. Faraday Soc. 33 (1937), 1272-1294. 

[11] S. Friedland, A proof of a generalized van der Waerden conjecture on per- 
manents, Lin. Multilin. Algebra 11 (1982), 107-120. 

[12] S. Friedland, On the entropy of Z-d subshifts of finite type. Linear Algebra 
Appl. 252 (1997), 199-220. 

[13] S. Friedland, Multi-dimensional capacity, pressure and Hausdorfi' dimen- 
sion, in Mathematical System Theory in Biology, Communication, Compu- 
tation and Finance, edited by J. Rosenthal and D. Gilliam, IMA Vol. Ser. 
134, Springer, New York, 2003, 183-222. 

[14] O.J. Heilman and E.H. Lieb, Theory of monomer-dimer systems, Comm. 
Math. Phys. 25 (1972), 190-232; Errata 27 (1972), 166. 

[15] D.S. Gaunt, Exact series-expansion study of the monomer-dimer problem, 
Phys. Rev. 179 (1969) 174-186. 

[16] J.M. Hammersley, Existence theorems and Monte Carlo methods for the 
monomer-dimer problem, in Reseach papers in statistics: Festschrift for J. 
Neyman, edited by F.N. David, Wiley, London, 1966, 125-146. 



34 



[17] J.M. Hammersley, An improved lower bound for the multidimesional dimer 
problem, Proc. Camb. Phil. Soc. 64 (1966), 455-463. 

[18] J.M. Hammersley, Calculations of lattice statistics, in Proc. Comput. 
Physics Con.., London: Inst, of Phys. & Phys. Soc, 1970, 1-8. 

[19] J. Hammersley and V. Menon, A lower bound for the monomer-dimer prob- 
lem, J. Inst. Math. Applic. 6 (1970), 341-364. 

[20] L.P. Hurd, J. Kari, and K. Culik, The topological entropy of cellular au- 
tomata is uncomputable, Ergodic Theory Dynam. Systems, 12 (1992), 255- 
265. 

[21] P.W. Kasteleyn, The statistics of dimers on a lattice, Physica 27 (1961), 
1209-1225. 

[22] C. Kenyon, D.Randall and A. Sinclair, Approximating the number of 
monomer-dimer coverings of a lattice, J. Stat. Phys. 83 (1996), 637-659. 

[23] M. Jerrum, Two-dimensional monomer-dimer systems are computationally 
intractible, J. Stat. Phys. 48 (1987), 121-134. 

[24] E.H. Lieb, The solution of the dimer problem by the transfer matrix 
method, J. Math. Phys. 8 (1967), 2339-2341. 

[25] E.H. Lieb, Residual entropy of square ice, Phys. Review 162 (1967), 162- 
172. 

[26] P.H. Lundow, Compression of transfer matrices. Discrete Math. 231 (2001), 
321-329. 

[27] J.F. Nagle, New series expansion method for the dimer problem, Phys. Rev. 
152 (1966), 190-197. 

[28] Z. Nagy and K. Zeger, Capacity bounds for the 3-dimensional (0, 1) run 
length limited channel, IEEE Trans. Info. Theory 46 (2000), 1030-1033. 

[29] L. Pauling, J. Amer. Chem. Soc. 57 (1935), 2680-. 

[30] K. Schmidt, Algebraic Ideas in Ergodic Theory, Amer. Math. Soc, 1990. 

[31] A. Schrijver, Counting 1-factors in regular bipartite graphs, J. Comb. The- 
ory B 72 (1998), 122-135. 

[32] I.M. Wanless, A lower bound on the maximum permanent in A^, Linear 
Algebra Appl. 373 (2003), 153-167. 

[33] W. Weeks and R.E. Blahut, The capacity and coding gain of certain 
checkerboard codes, IEEE Trans. Info. Theory 44 (1998), 1193-1203. 



35 



